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Abstract 

This work presents a detailed characterisation of the MAST Scrape Off Layer in L-mode. Scans 
in line averaged density, plasma current and toroidal magnetic field were performed. A compre¬ 
hensive and integrated study of the SOL was allowed by the use of a wide range of diagnostics. In 
agreement with previous results, an increase of the line averaged density induced a broadening of 
the midplane density profile. This increase was not correlated with divertor detachment, as con¬ 
firmed by the systematic increase of the target ion flux and decrease of the D^/Dq, emission. Also, 
no clear correlation is found with the density of the neutral particles at the wall. At comparable 
density levels, discharges with higher current did not show broadening. Outer target ion saturation 
current and heat flux decay lengths were measured and compared with midplane data. For the 
saturation current, the upstream projections of the target values, based on diffusive models, did 
not match the midplane measurements, neither in amplitude nor in trend, while agreement was 
found for the heat fluxes, suggesting a different perpendicular transport mechanism for the two 
channels. Furthermore, the value of the heat flux decay length was quite insensitive to changes 
in the thermodynamic conditions, in agreement with recent scaling laws. In all the cases studied, 
sawtooth oscillations were present but they simply rescaled self-similarly the target profiles. The 
separatrix conditions changed significantly during a sawtooth cycle, but the heat flux decay length 
and divertor spreading factor remained nearly constant, indicating that these quantities are rather 
insensitive to the upstream thermodynamic state of the SOL. 
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I. INTRODUCTION 


As magnetic fusion research progresses towards reactor relevant conditions, it becomes 
clearer and clearer that the plasma exhaust and its consequent surface interaction will 
strongly constrain the o^rational space and determine whether the next generation ma¬ 
chines will be successful [l|. The level of interaction of the plasma with the solid structures 
surrounding it is largely determined by the properties of the Scrape-Off Layer (SOL), a nar¬ 
row region outside the separatrix which is magnetically connected with the divertor target 
or the walls of the machine j^. One of the defining features of the SOL is its width, often 
described as the decay lengths of some combination of thermodynamic quantities. 

Of particular interest are the particle decay length, = {d\ogn/dr)~^, which captures 
the intensity of the particle fluxes towards the walls and Xq = {d\ogq\\/dr)~^, which de¬ 
termines the divertor area over which the power is released (here n is the plasma density, 
q\\ the parallel heat flux and r the radial coordinate). Due to the large difference between 
parallel and perpendicular transport, the latter being much less efficient, the SOL width is 
significantly smaller than the machine size, i.e. Xq ^ R and A„ R, where R is the major 
radius. 

Despite the small size of the SOL, the particle flux at the wall, mediated by filamentary 


structures 3|, 1^, can still be signihcant. As several experimen 
years, this is particularly true in high density operations js, 6 


rave shown over the last 15 


in| when a shoulder forms 


in the density profile. Importantly, if ions in the filaments maintain an energy above the 
guttering threshold when they reach the solid surfaces they can induce localised erosion 
thus degrading the machine and releasing impurities that pollute the plasma. Also, 
the density prohle in the SOL affects the coupling of the RF antennas, thus determining 
the efficiency of the heating deposition and it affects the plasma fuelling by influencing the 
neutral penetration in the core. 

Because of the small size of the SOL, the surface interaction area at the divertor is 
comparably small, thus leading to unacceptable transient and steady state heat loads in 
high power machines, which entails the need for control techniques (e.g. tailoring of the 
divertor geometry, impurity seeding, advanced conhgurations). Similarly, volumetric losses 
(such as radiation) which would alleviate the divertor loads are limited by the SOL volume 
which is small compared to the total plasma volume. 
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In the last few years, an important experimental effort was devoted towards the empirical 
understanding of the scaling of the heat flux decay length with respect to plasma parameters. 
This work, based on infra-red thermograpy, resulted in multi-machine regression analyses 
which suggested that in attached H-mode Xg (at the outer midplane, extrapolated from 
the divertor data) mainly depends on the inverse of the poloidal magnetic held lll-ll3|. 
Interestingly, following L-mode studies showed remarkably similar trends. 

It is important to note that the SOL width and its features depend on the position, as 
one moves along the held lines from upstream (i.e. at the midplane) to downstream (i.e. at 
the target) 1^. The decay lengths change because of the magnetic conhguration (i.e. the 
hux expansion) and because of the variations in the cross held and parallel transport, which 
depend on the plasma parameters. 

In this paper, we report a detailed characterisation of the L-mode Scrape Oh Layer in 
the Mega Ampere Spherical Tokamak (MAST). Our study benehts from a wide range of 
diagnostics, which allowed us to correlate the outer midplane SOL to the divertor region, 
thus providing useful insight in the parallel ehects. In addition, we extended our investiga¬ 
tion to the inner divertor region, which behaves in a remarkably diherent way with respect 
to its outer counterpart. We devoted particular attention to the investigation of the mech¬ 
anisms that induce the density prohle broadening upstream and how this translates in the 
downstream region. 


II. EXPERIMENTAL CONDITIONS AND DIAGNOSTICS 


All the results presented here were obtained in MAST 1^, which is a small aspect 
ratio machine with major radius R ^ 0.9m and minor radius a ~ 0.6m. Other important 
peculiarities of the machine are the particularly open divertor and the large distance between 
the separatrix and the low field side wall (~ 0.5m). It is also useful to remember that the 
MAST solenoid fringing field causes the outer strike point to sweep the divertor target by 
tens of centimetres during the plasma discharge. 

We investigated Ohmic L-mode plasmas in attached conditions over a wide range of 
upstream SOL conditions. All the discharges were performed with double null diverted 
configurations. It is important to remark that, due to the particular divertor and wall 
geometry of MAST, detachment is particularly difficult to achieve 0 and none of the 
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FIG. 1. Line averaged density and ^95 for the reference [(a) and (b)], high current [(c) and (d)]and 
low magnetic field [(e) and (f)]cases. 


results presented here was in that regime. 

Our reference case has Bt = 0.585T and Ip = 400fcA, while the line averaged density 
during the flat top was varied between three levels, Ue ~ [1-3,1.65, 2.2] x This was 

done by varying the deuterium injection from piezo valves situated in the lower and upper 
divertor regions and at the high held side wall. We also performed a set of density scans at 
higher current. Ip = 600kA, with reference held and one at low magnetic held, Bt = 0.4T, 
with reference current, for a total of eight diherent experimental conditions (the low density, 
low magnetic held case was not suitable for the analysis due to large MHD instabilities). In 
FiglUthe line averaged density (measured by the MAST interferometer) and the edge safety 
factor, g 95 , are shown for the reference, high current and low hied cases. 

The divertor region is well monitored with 570 hxed Langmuir probes (FLP) embedded 
in the centre column (where the inboard strike point lies) and target structures. The probes 
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cover the upper and lower strike points and are 3mm apart on the inboard side and 10mm 
on the outboard (the time resolution reaches a maximum of IkHz). In addition, the target 

n 

heat fluxes are measured with two Infra Red cameras (IR) [18| . The hrst looks at the upper 
strike points (inboard and outboard) in the long wave infra red (7.6-8.9/im) with a spatial 
and temporal resolution of 6.26mm and 2.bkHz, while the second looks at the lower strike 
points (inboard and outboard) in the medium wave infra red (4.5-5pm) with a spatial and 
temporal resolution of 5.9mm and 0.33kHz. 

Upstream, electron density and temperature are measured by a high resolution Thomson 
scattering system (HRTS) [l9j capable of resolving the core and edge profiles at 0.2AkHz 
with 10mm radial resolution. Furthermore, the ion saturation current and electrostatic 
plasma potential were monitored by a mid-plane reciprocating probe with a Mach head 


mounted on it 


20|. During each discharge, only one probe reciprocation was performed. 


In addition, spectrometers measured the Da (656.Inm), (434.Inm) and Cu (514?7,m) 
light emission integrated over tangential and radial views at the outer mid-plane and in the 
divertor regions. The sampling frequency of these measurements was 50kHz. To complement 
these measurements, a Da linear camera, which is a 1024 element 2D CCD array, provided 
the radial profile of the emission at the midplane. Finally, we used a standard Bayard- 
Alpert ion gauge, absolutely calibrated and equipped with a fast data acquisition system, 
to measure with a time resolution of 0.1 ms the pressure of neutral deuterium at the outer 
midplane. It should be noted that, since the ion gauge is located in a recess shielded from 
the tokamak magnetic field, its response has a delay of a few milliseconds with respect to 
the variation of the neutral pressure in the vessel, due to the finite conductivity of the pipe. 


III. PHENOMENOLOGICAL DESCRIPTION 

In this Section, we discuss the general features of the eight cases analysed, in order 
to better interpret the quantitative data associated to the SOL decay lengths. All the 
measurements discussed below were taken in the time window 0.25 —0.3sec unless otherwise 
stated. In this period, for each discharge the main plasma parameters had only small 
variations, e.g. the relative change between highest and lowest value was between 10 — 20% 
for g 95 and between 5 — 10% for Ue. The position of the separatrix at the midplane, where the 
upstream temperatures are evaluated, is estimated using an optically constrained magnetic 
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reconstruction obtained with the EFIT code 


2l| . which typically has a Icm accuracy. 


For convenience, the discharges are clustered in ’’reference”, ’’high current” and ’’low 
held” sets. In the last two subsections we compare the discharges in order to provide a more 
integrated overview of the plasma conditions, in particular at the divertor and low held side 

wall. 


A. Reference set 


At the low density level, the Greenwald fraction is few = 0.36. The upstream electron 
temperature obtained from the HRTS is ~ 20eV, as conhrmed by more precise mea¬ 
surements obtained with a retarding held energy analyser in a similar discharge 22| , which 
also suggest a separatrix ion temperature twice as high. The weak thermal coupling between 
ions and electrons is due to the low collisionality of the upstream SOL. The power crossing 
the separatrix, Psol, obtained by subtracting the time variation of the energy and the core 
radiated power measured by the MAST bolometer from the Ohmic power, is between 350kW 
and 400kW during the hat top. Using a simple two point model to estimate the upstream 


temperature, we obtain: Tg 


17/2 


PsolL\\ 


T -I- 

4:2TrR\q{Bp/BT)K.0,e 


20eU, in agreement with 
our direct measurements. In this calculation we have used {Bp/Bt) ~ 0.3 and Ri 0.02m 


2/7 


as discussed in Section IV Bl and is the temperature independent part of the parallel 
electron conductivity. Also, this temperature agrees very well with the empirical scaling 


Tg^j = 1.4 X (see Eq.4.42 of j^), and so do almost all the other seperatrix 

temperatures given in this and the following Subsections. Since the HRTS data show a 
separatrix density ~ 0.35 x 10^®m“^, and the connection length, L||, is around 13m, the 
electron collisionality can be estimated as i 2 ^, g = 10“^®ng^„L||/Tg^„ Ri 11, which would imply 
a sheath limited regime, but close to the transition to the conduction limited regime. In the 
outer upper divertor, the FLP measured a peak density Ue^t ~ 0.35 x 10^®m“^ and an electron 


temperature at the strike point position Tg j Ri 14eU. Based on the results of 22|, l26j, we ex¬ 
pect an equilibration between target electron and ion temperature, which should be roughly 
the same, i.e. Ti^t ~ Te^t- Comparing upstream and downstream data, we hnd pressure con¬ 
servation along the held lines, which can be formulated as: rie^uiTe^u + Ti^u) ~ 2ng^i(Tg_t-|-Tj i) 


(see e.g. 


2d|). 


The second density level corresponds to few = 0.48. The higher line averaged density 


6 


















corresponds to an increase in the separatrix value which is ne,u ~ 0.49 x and 

also leads to a cooling of the upstream temperature, which is around ~ lAeV. The 
conditions at the upper outer target are rze,* ~ 0.7 x and ~ 7eV, thus showing 

again pressure conservation in the assumption of similar ion/electron temperature ratios as 
in the previous case. The upstream electron collisionality is z/* e ~ 32, which would suggest 
a conduction limited regime, compatible with the fact that the target temperature is half 
the upstream. 

The third level, at few = 0.6, shows a similar upstream temperature, ~ lOeld, as 
the intermediate density case, which is reasonable considering a conduction limited regime 
and comparable power crossing the separatrix, while the separatrix density reaches ~ 
0.55 X This is again compatible with the measurements in 22[, which also indicate 

an upstream ion temperature a factor two higher. At the target, Ue^t ~ 0.75 x is 

measured, while the temperature data show a large scatter. Based on the results in |22| , we 
estimate the electron and ion target temperature for this case at around 5eV. Like in the 
previous cases, pressure is conserved and the upstream collisionality for this density level is 




71. 


B. High current set 

The second set of discharges was performed at a higher current {Ip ~ GOOkA) and at the 
reference magnetic held. The main consequences of this change are a larger Ohmic heating, 
Pohm ~ Ip, which leads to a temperature increase with respect to similar density levels in 
the reference cases and a reduction of the edge safety factor, qg^ ~ (this also shortens 
the connection length to roughly 85% of the reference value). Figures [T] and [2] show that 
these effects are observed while the core density prohles are reasonably aligned to each other 
apart from the low density level, in which the line averaged density is around 20% lower 
than the reference discharge. 

In the low density level, with few ~ 0.19, the measured electron density and temperature 
at the separatrix are Ue^u ~ 0.28x10^® and ~ 22eV. Due to the low collsionality, ~ 6, 
the divertor is expected to be in a sheath limited regime, which is conhrmed by the fact 
that target electron temperature is similar to the upstream value 20eV. For this case 

we do not have ion temperature measurements but we expect low thermal coupling between 
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FIG. 2. Time averaged profiles of the electron temperature (upper row) and density (lower row) 
for the 8 discharges during the flat top. Only the outboard half is shown. From left to right, the 
columns show increasing density levels. Blue, red and green curves represent the reference, high 
current and low held sets. In the lower row, the horizontal dash-dot lines represent the values of 
the line averaged density of the different discharges. 


ions and electrons, due to the low collisonality. Given the measured electron density at the 
target, rig,* ~ 0.3 x 10^®m“^, pressure conservation, which should be a consequence of the 
sheath limited regime, would be satisfied with Tj rs 3Te_„ if Ri Tg^i, which is comparable 
to the values of the reference case. 

The second density level, at few ~ 0.32, is characterised by ~ 0.51 x 10^® and 
~ 21eV, which leads to an upstream collisionality z/* g Ri 9. The target data give 
rig t Ri 0.7 X 10^® and Tg_t Ri 15eV. Pressure conservation would again be satished if the 
upstream ion temperature was roughly three times higher than the electron temperature 
under the condition that the target thermal coupling was strong. 
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Finally, the last density level corresponds to few ~ 0.41. Upstream and downstream 
measurements give rie^u ~ 0.6 x 10^®, ~ 12eU, ne,t ~ 1-3 x 10^® and Te^t ~ 7eU. With 

these conditions, the divertor should still be in attached conditions and conduction limited, 
since z/* e ~ 45. Pressure conservation would require an upstream ion temperature around 
twice as high as the electron, compatible with higher collisionality with respect to the lower 
density levels and hence higher thermal coupling. 


C. Low field set 

The final set of discharges, which does not include a low density case, was performed at 
a lower magnetic field, Bt = 0.4T, and at the reference current. This gave a ggs similar 
to the high current cases since ggs ~ Bx/Ip and an Ohmic heating similar to the reference 
discharges. The core density profiles match reasonably well the reference and high current 
sets, while the temperature profiles are comparable to the reference cases, albeit lower in 
the innermost part of the plasma, see Fig 121 

As expected, the Greenwald fraction of the intermediate density discharge, few = 0.44, is 
similar to its reference case counterpart, and so are the separatrix density, ~ 0.5 x 10^®, 
and temperature, ~ 12eU. The associated collisionality, which takes into account 
the reduction of the connection length due to the lower ggs, is z/*_e ~ 24. At the target, 
~ 0.6 X 10^® and ~ 9eU. 

Also in the high density case the matching with the reference case is good as upstream we 
measure few = 0.63, ne,u ~ 0.7 x 10^® and Te^u ~ lOeU, which gives z/*^e ~ 49. Downstream, 
the Langmuir probes give rie^t ~ 0.7 x 10^® and ~ 7eU. 

D. Radiation, total ion flux and divertor regime 

Using the combined data of the discharges, we can shed some light on their divertor 
regime and in particular on whether detachment occurred. The first indication of attached 
conditions comes from the fact that pressure balance seems to be satisfied in all the dis¬ 
charges (see Subsections above), which indicates that only a small amount of momentum 
is transferred to neutral particles. This estimate is, however, very rough as measurements 
are sometimes incomplete and subject to significant error bars. While pressure conservation 
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should be carefully checked in the future, we present in the following more solid arguments 
which suggest absence of detachment. 

For each set of discharges (reference, high current and low held), during the hat top phase, 
a higher line averaged density corresponds to higher core and edge radiation, measured by 
the bolometer (not shown) and the deuterium and carbon line emission, see FigJHl At 
the same time, the D^/Da ratio at the target steadily decreases, as expected when the 
divertor density increases in conditions in which excitation dominates over recombination, 
i.e. before detachment 23(]. This can be clearly seen in the plots on the second row of Figj3] 
when diherent density levels of the same set of discharges are compared. It is worthwhile 
remarking that a reversal of this trend in the reference set was transiently observed in a very 
high density shot, corresponding to an average fie ~ 2.9 x This discharge could not 

be included in the present analysis as it showed large oscillations in the plasma parameters 
and no steady state period due to the difficulty of feedback controlling the density. 

Coming back to the analysed discharges, the total ion flux to the target, proportional to 
the radial integral of ) ~ Jsat,t-, was estimated using the divertor Langmuir probes. 
Larger ion fluxes were observed when the line averaged density was higher, again confirming 
that the discharges were not detached. This is shown in FigJU where the total ion flux, 
Fj ~ J Jsat,tRdR is plotted as a function of the upper outer strike point position during 
the analysed time window. The reduction of the total ion flux as RstrikePoint increases is 
attributed to the lower plasma temperature at the target as the divertor leg moves outward. 


E. Neutral density and plasma sources 

The neutral pressure at the wall, measured near the midplane by the fast ion gauge, 
allows us to estimate the mean free path of the neutrals and their likelihood of being ionized 
before reaching the separatrix. In the discharges presented here, the neutral pressure showed 
a good correlation with the line averaged density and appeared to be quite insensitive to 
other operational conditions, e.g. the plasma current or the toroidal magnetic held, see the 
lower panel of FigJS] 

More details on the plasma sources and the neutral density can be obtained by comparing 
the radially integrated Da emission at the midplane and the neutral density at the wall. In 
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FIG. 3. Radiation measurements for the 8 discharges. The first row shows the emission, the 
second the D^/Da ratio, the third the Cjj emission. From left to right, the columns show the 
reference, high current and low magnetic field discharges. The signals were smoothed with a box¬ 
car average with a 10ms window to filter out the sawtooth oscillations. All the measurements were 
taken from radial views at the upper divertor. 


all the cases treated, the emission is localized around the outer separatrix, as the linear 
camera measurements show. The radially integrated emission measurements (see FigIS]) are 
therefore representative of that region of the plasma. As FigIS] shows, for both separatrix 
and wall signal, the strongest dependence is on the plasma line averaged density, suggesting 
that the particle sources are similar between sets at comparable Ug. 


It is also useful to mention that previous work 2^ has shown that, for a typical discharge 
duration, the MAST wall is expected to retain almost all the particles impinging into it. 
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FIG. 4. Integrated ion saturated current measured by the Langmuir probes, used as a proxy for 
the total particle flux at the outer target as a function of the strike point position during the flat 
top phase. Different symbols represent different discharges. 


IV. UPSTREAM PROFILES AND DECAY LENGTHS 

A. Ion saturation current 

The plasma density and the ion saturation current are strongly related. A Mach probe 

mounted on a reciprocating arm was used to evaluate the upstream prohles of the latter. High 

Resolution Thomson Scattering data were in good agreement with the RP, as shown in Figl6] 

1 /2 

for the reference discharges, with the ion saturation current estimated as Jsat,u ~ ne,uTe,u . 

We remark that in the evaluation of the prohles, we processed only data during the hrst 
part of the reciprocation, i.e. before the probe started to move out of the plasma. This was 
necessary as the probe reciprocation perturbs the local properties of the plasma so that on 
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FIG. 5. Upper panel; time evolution of the midplane radial Da signal smoothed with a moving 
average technique (using a 10msec window). Lower panel: time evolution of the neutral density 
at the wall. The horizontal dotted lines show the separation between the three density levels. The 
oscillations in the signals are attributed to the response of the plasma density control system. 


its way out, the probe measures steeper profiles. Importantly, if both inward and outward 
data were averaged and processed together, the flattening in the far SOL would disappear. 
In practical terms, our data were collected between 0.23 —0.28sec but with different starting 
points and window length depending on the discharge (all the time windows had a minimum 
length of 30msec). 

The comparison between the different discharges shows that the SOL prohle broadening 


is consistently correlatec 


in other machines 


m, 


with the increase of the line averaged density, as already observed 


10| . see FigITl However, there is a qualitative and quantitative 


difference between the reference set and the other discharges. In the former, the slope of 
the prohles is steeper close to the separatrix, while it flattens in the far SOL, especially for 
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FIG. 6. Comparison between the normalized ion saturation current HRTS averaged profiles (dashed 
lines) and RP profiles (thick line). A thin line represents the exponential fit to the RP data used 
to calculate Xj^u- The points indicate the position and value of the HRTS channels, with an empty 
symbol representing less reliable channels that had more than 50% unusable data in the period of 
the average. The gray area (plotted only for the reliable channels) represents the uncertainty in 
the HRTS measurements. 


large Ug. Also the high density low magnetic field case showed a similar trend. There is 
extensive experimental literature documenting this two region feature of the diverted SOL 


25 


28| . although the mechanisms behind this phenomenon remain elusive. In the high 


current cases, however, the profiles did not show any visible flattening in the far SOL. This 


is in agreement with the results reported in 


25| for a closed divertor, which suggests that the 


current effect is insensitive from divertor geometry and closure. In addition, in most of the 
cases analysed, the gradients of the profiles in the near SOL and just inside the separatrix 
did not show sharp discontinuities. 
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FIG. 7. Profiles of the ion saturation current normalized to its separatrix value as a function of 
the distance from the separatrix. The profiles from the RP measurements are smoothed using a 
moving average with a bin of 15msec. For each set, the symbols (v), (o) and (a) represent the 
low, intermediate and high density levels respectively. Straight lines, representing the exponential 
profiles with decay lengths given in Table [H are plotted for comparison. 


These reciprocating probe results were confirmed by measurements taken with the HRTS 
system, which also reiterated that the density follows the ion saturation current profile 
[HI The data shown were obtained by accumulating the scattered signal time traces in a 
stationary plasma over a period of 50ms. Since the lasers fire at 240Hz, this corresponds to 
12 scattered pulses. Following the accumulation process, the signal time traces are htted to 
obtain the number photons in each spectral bin of a polychromator, which are in turn fitted 
to obtain density and electron temperature. While the diagnostic typically measures well 
down to 2 — 5 X in single pulse mode, accumulating over 12 laser pulses allows us to 

measure with conhdence down to ~ 0.6 — 1.4 x 10^®m“^. It is important to accumulate the 
time traces instead of htting each scattered signal individually and averaging the integrals, 
as the signals for an individual time trace are close to the noise floor. It should be noted 
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FIG. 8. Profiles of the density and ion saturation current normalized to their separatrix value as a 


function of the distance from the separatrix. The prohles are obtained from HRTS measurements 


with the technique described in the text. 


TABLE I. Ion saturation current decay lengths calculated upstream in the near SOL (in cm). 



low He 

intermediate Ug 

high rie 

Reference 

3.58 ± 0.01 

5.72 ± 0.01 

7.41 ± 0.04 

High Current 

4.05 ± 0.01 

/ 

7.82 ± 0.02 

Low field 

/ 

2.58 ± 0.03 

6.27 ± 0.01 


that our diagnostic has a minimum measurable temperature of approximately 5eV due to 
the proximity of the final spectral bin to the laser wavelength. 

From these measurements, we observe that the decay length in the near SOL ranges 
between ~ 3.5cm and ~ 7cm in all the experimental conditions examined. The values of Aj 
with their fit errors are given in Tabled! The absence from the table and from FigJTjof the 
intermediate density in the high current set is due to the fact that the reciprocation was not 
deep enough to evaluate the near SOL features. 
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From these results, the dominant parameter that controls the near SOL ion saturation 
current decay length seems to be the line averaged density. The comparison between the 
reference and the high current sets also suggests low sensitivity to Jp, and hence that the 
Greenwald fraction is not a good scaling quantity for Aj. 

In the far SOL, Aj responds to variations of rTg more rapidly than in the near SOL for 
the reference discharges. Also, different^ from the near SOL, plasma current affects the 


gradients in this region. Previous work [5 


suggested that an increase of the ionization 


sources due to the neutral recycling at the wall could provide a possible explanation for 
the far SOL broadening. However, the results of Sec lIII El suggest that the level of neutrals 
depends exclusively on the line averaged density, while the broadening has a clear dependence 
on Ip and Bt as well. This implies that wall recycling cannot explain the flatter prohles in 
the far SOL of the reference set. 

In our experiment, we do not observe a connection between broadening and detachment, 
which is sometimes believed to affect the SOL prohles. While the divertor plasma is changing 
between different discharges, so is the SOL above the X-point, thus making it extremely 
difficult to pinpoint the mechanism for the prohle battening to a particular SOL location, 
which might be due to global properties of the edge plasma. 


B. Parallel heat flux 

In our experiments, no direct measurements of the heat flux profile are available upstream. 
However, the decay length close to the separatrix can be estimated by extrapolating to the 
SOL region reasonable combinations of the HRTS and RP data. This is justified by the 
continuous nature of the gradients across the separatrix observed in the ion saturation 
current. As a standard approximation j^, if conduction is the dominant parallel heat loss 
mechanism, the heat flux decay length is estimated as a fraction of the electron temperature 
decay length, A^ ~ 2/7Xt- On the other hand, if parallel convection is assumed to be more 
important, a better approximation is A“^ ~ A“^ + 3/2Ay^ (from gy ~ Both regimes 

constitute a simplification of the actual parallel heat loss physics, which at the very least is a 
combination of both effects. In principle, in perfectly symmetric double null configurations, 
the parallel heat flux at the midplane should be exactly zero and therefore Xq is an ill 
defined quantity. Therefore, our measurements are simply reflecting combinations of the 
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TABLE II. Parallel heat flux decay lengths calculated in the near SOL (in cm). Values extrapolated 
from HRTS data. 



low rie 

intermediate ng 

high rif. 

Reference (cond.) 

0.62 ± 0.12 

0.90 ± 0.49 

1.60 ± 0.28 

Reference (conv.) 

1.35 ± 0.28 

2.03 ± 1.20 

3.19 ± 0.61 

High Current (cond.) 

0.85 ± 0.19 

/ 

1.17 ± 0.59 

High Current (conv.) 

1.72 ± 0.42 

/ 

2.69 ± 1.45 

Low field (cond.) 

/ 

0.95 ± 0.17 

1.50 ± 0.55 

Low field (conv.) 

/ 

1.45 ± 0.29 

2.86 ± 1.15 


temperature and density profiles and as such they should be interpreted. However, it makes 
sense to use the convective approximation to compare the upstream with the downstream 
data, since the sheath at the target enforces q\\^t ~ |^. In other words, the convective 

estimate for our upstream measurements is representative of the one direction heat flux that 
would be obtained if a solid surface was present at the midplane. 


For the sake of completeness, we give here estimates of Xq in both limits, knowing that 
these values are only indicative, although probably in the right order of magnitude, and 
that Xq convective is probably more relevant. To improve the accuracy of the estimate, 
instead of using A„ in the calculation of the convective decay length, we use Aj, so that 
A“^ Ri Xj^ + A^h The upstream decay lengths calculated with this method are summarized 
in Tabim The errors in the table are a result of the fitting procedure (based on the last 5 
channels of the HRTS) and were calculated assuming a 95% conhdence interval. Note also 
that the temperature decay lengths from the HRTS can be easily extrapolated from the 
conductive Xq. 


The upstream near SOL heat flux profile appears to have a modest dependence on the 
plasma conditions, with Xq around 1.5-3cm for the convective model and 0.5-1.5cm for the 
conductive model. This is compatible with the observation that the dominant dependence 
of Xq i s on the magnetic configuration, and in particular on the poloidal magnetic field 


11 


, ll3|, [l^. However, also the line averaged density, or possibly the collisionality, seem to 


have a moderate but visible effect on Xq, which somewhat increases with rig or 
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V. TARGET MEASUREMENTS 


The radial profiles of both the target ion saturation current and parallel heat flux were 
measured directly with Langmuir probes and infra red thermography respectively. All the 
outer divertor data were analysed using the same procedure, based on a htti^ function that 
assumes a diffusive broadening of exponentially decaying upstream profiles jlll. l33l | . Its form 
is given by: 


/(^) 



X — Xq 


+ fb, 


'gi 


( 1 ) 


where / can represent the parallel heat flux or the ion saturation current, Sf is the divertor 
spreading factor and Fexp is the target flux expansion. The htting parameters /o, fbg and xq 
measure the amplitude of the field, its background value and the position of the strike point. 
Note that, in this model. A/ represents the projected upstream decay length, which we use 
as a reference throughout the paper and is directly comparable with the data presented in 
Sec lIVl To assure compatibility between upstream and downstream measurements, we have 
restricted the analysis of the target data to the same time window defined in Sec lIVl even 
though the divertor diagnostics were active throughout most of the discharge time. 

Figure [9](a) shows a typical prohle of the parallel heat flux (averaged over transien 
caused by filaments) and its fitting curve which, as also observed in other machines 
usually provides a very good match. This is not the case, however, for the prohles in the 
inner target, which require a different approach, as discussed in Sec IV Cl 


: events 
12, 131, 


A. Ion saturation current 

To increase accuracy, the analysis of the Jsat,t profiles was restricted to measurements 
with a relative error below 15%. The natural sweep of the strike point was used to construct 
a better resolved profile by overlapping the profiles at different time frames, after correcting 
for their increasing radial shift. This was done during the flat top phase and in a 30-50 
msec window. The upper outer divertor decay lengths and divertor spreading parameters 
obtained with this method are summarized in Tables m and lYl 
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(a) 




FIG. 9. Heat flux profiles at the upper outer (a) and inner (b) targets for discharge 28996 at 
t = 0.2508 sec. The markers are the IR measurements, while the curves are fits in the form of 
Eq(T]for (a) and of Eql2]for (b). R is the major radius at the divertor plate and Z is the vertical 
coordinate on the central solenoid. Dotted lines represent the position of the strike point. 


B. Parallel heat flux 


The good temporal and spatial resolution of the IR system allowed us to treat these data 
in a different way with respect to the ion saturation current. In particular, we htted each 
time frame, thus obtaining the detailed time evolution of both Ag and Sq. 

It is important to remark that in all our discharges we observed small sawteeth oscilla¬ 
tions. This has interesting consequences for the deposition pattern at the divertor targets, 
producing a dynamic behaviour in several physical quantities, including Xg. A detailed dis¬ 
cussion of these effects will be presented in a companion paper [3^. Despite the transient 
behaviour induced by the sawteeth, it was possible to extract measurements representative 
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TABLE III. Ion saturation current decay lengths, Aj, calculated at the upper outer divertor (in 
cm). Fitting errors calculated using a 95% confidence value. 



low He 

intermediate Ue 

high rie 

Reference 

1.37 ± 0.03 

2.71 ± 0.03 

2.77 ± 0.05 

High Current 

1.54 ± 0.04 

1.79 ± 0.03 

2.57 ± 0.06 

Low field 

/ 

1.46 ± 0.02 

2.62 ± 0.02 


TABLE IV. Same as Tab lllll for the ion saturation current divertor spreading parameter, Sj (in 



low He 

intermediate Ue 

high rie 

Reference 

0.74 ± 0.02 

0.53 ± 0.01 

0.62 ± 0.02 

High Current 

0.88 ± 0.03 

0.52 ± 0.02 

0.61 ± 0.02 

Low field 

/ 

0.66 ± 0.01 

0.56 ± 0.01 


of steady state condition. Indeed, while the upstream density and especially temperature 
show an oscillating behaviour, the deposition prohles remained self-similar during the saw¬ 
tooth cycle, apart from a relatively short period of time during the crash (~ 1msec). This 
suggests that both Xg and Sg have a low sensitivity to the thermodynamic properties of the 
upstream plasma. 

This is conhrmed by comparing the decay lengths of the different discharges performed, 
all of which show a similar value of Xg. For comparison, the average Xg and Sg, obtained 
by removing the data corresponding to the crash phase of the sawtooth, are summarized in 

Tables El and El 

TABLE V. Parallel heat flux decay lengths, Xg, calculated at the upper outer divertor (in cm). 
Fitting errors calculated using a 95% confidence value. 



low rie 

intermediate Ue 

high rie 

Reference 

1.83 ± 0.12 

1.90 ± 0.17 

1.96 ± 0.22 

High Current 

1.88 ± 0.16 

1.97 ± 0.08 

2.13 ± 0.15 

Low field 

/ 

2.10 ± 0.12 

2.02 ± 0.20 
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TABLE VI. Same as TablVlfor the heat flux divertor spreading parameter, Sq (in cm). 



low He 

intermediate Ue 

high rie 

Reference 

0.32 ± 0.04 

0.81 ± 0.08 

0.76 ± 0.11 

High Current 

0.36 ± 0.06 

0.35 ± 0.03 

0.70 ± 0.07 

Low field 

/ 

0.26 ± 0.03 

0.92 ± 0.10 


TABLE VII. Parallel heat flux divertor spreading parameters calculated at the upper inner divertor 
(in cm). Eitting errors calculated using a 95% confidence value. 



low He 

intermediate He 

high rie 

Reference 

0.98 ± 0.09 

0.79 ± 0.08 

1.28 ± 0.20 

High Current 

1.03 ± 0.08 

0.99 ± 0.06 

1.02 ± 0.05 

Low field 

/ 

0.85 ± 0.06 

0.76 ± 0.12 


C. Inner divertor measurements 

In the double null configurations that we have analysed, only roughly 10% of the power 
crosses the separatrix towards the high field side. In the low power environment of the inner 
target, the IR prohles show a symmetric conhguration with respect to the strike point, see 
Figj9](b). A consequence of this is that when the Wagner/Eich function, Eq{Tl is fitted to the 
inner target IR data the decay length is poorly constrained (it would lead to large errorbars). 
This suggests that, differently from the outer target, Xq is smaller than Sq, and hence the 
prohle is mostly determined by the mechanisms below the X-point. In our analysis, we fitted 
the prohles with a simple Gaussian function of the form: 

/(x)=/„e’('V‘) (2) 

which provided a better £t to the data. Table IVIII summarizes the measurements taken. 


VI. COMPARISON BETWEEN UPSTREAM AND TARGET MEASUREMENTS 

Comparing upstream and downstream measurements allows us to gain some insight into 
the mechanisms that transfer particles and energy towards the divertor. In particular, the 
downstream data are now commonly interpreted with analysis techniques that assume cross 
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field diffusion 


33| from a midplane profile (interestingly, they also postulate a convective 


parallel transport for the energy). Figure [10] shows that midplane and target measurements 
are globally in good agreement for the heat flux decay length, while they are not for Aj. 

These results question the validity of the standard diffusive approach for the particle 
cross held transport in the SOL, which is likely to be governed by the complex motion 
of high density hlamentary structures that are observed in all our discharges. It is also 
interesting to notice that the near SOL Aj has a well dehned increasing trend with the 
upstream collisionality, both at the midplane and at the target. 

In Figiini the heat hux decay lengths upstream are evaluated with convective and conduc¬ 
tive assumptions, with the former displaying a better agreement with the projected target 
data for low collisionality and vice versa for the latter. Overall, both upstream estimates 
show an increase of Xq^u with the collisionality. On the other hand, the projected Xg is 
remarkably constant in different experimental conditions, which suggests that this quantity 
has a weak dependence on both thermodynamic and magnetic properties of the plasma. In 
addition, it seems like the upstream density broadening does not affect signihcantly the heat 
flux target prohles. A proper estimation of the scaling factors would require a database 
extended to more than 8 discharges (e.g. see jl3|). 

In this respect, it is useful to compare our measurements with recent experimental or 
theoretical scaling laws of Xq. The multi-machine regression based on infra red analysis 


performed in 


14| predicts a decay length of 1.4-1.9cm for our dataset, but with large error 


bars that can go one order of magnitnde above or below the nominal value. The theoretical 


scaling of 


30| . obtained with midplane 2D turbulence simulations, provides a good £t to 


the the convective npstream measurements, especially considering the low errorbars on the 
theoretical predictions (corresponding to a relative error of 15%). An exception is the low 


density high current case, which is characterised by the 
This might suggest that the flux limiter closure used in 


owest collisionality of the dataset. 


parallel dynamics and non-local approaches 


30| is not snfficient to captnre the 


3l| might be needed to improve the agreement. 


However, the matching of this midplane model is not as satisfactory where the projected 
target data are concerned. An heuristic theoretical model based on drifts was proposed 


m 


32| and its predictions give a Xg around 0.9-1.45cm. An important caveat is that the 


Goldston model j32|] was developed for H-mode plasmas. A snmmary of these comparisons 
is shown in FigJTTl 
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FIG. 10. Parallel heat flux (left panel) and ion saturation current (right panel) near SOL decay 
lengths as a function of the electron collisionality. The symbol (•) represents the upstream projec¬ 
tion of the target measurements, (o) the upstream measurements (with convective assumption for 
the heat flux) and (x) the upstream measurements with conductive assumption for the heat flux. 
The solid curves on the right panel show the linear fit of the upstream and target projected data. 

VII. SUMMARY AND CONCLUSIONS 

In this paper, we discussed a series of experiments carried out in the MAST spherical 
tokamak aimed at characterising the behaviour of the SOL decay lengths as a function of the 
plasma conditions. Our results show that the density and parallel heat flux profiles respond 
to changes in the plasma in a different way, both at the midplane and at the target. 

Only the particle exhaust, for example, is significantly affected by the density (or the 
fuelling), a higher level of which corresponds to broader density profiles in the SOL. Such 
broadening occurs both in the near SOL, where the exponential decay becomes less steep. 
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FIG. 11. Left panel: comparison between the parallel heat flux measured upstream and the scaling 
law of 30| (o, theoretical). Right panel: comparison between the upstream parallel heat flux 
projected from the target measurements and the scaling laws of [l^ (o, experimental) and of 321 
(x, theoretical), the solid lines are the errorbars on the measurements or on the predictions. The 
thin dashed lines represent the diagonals. 


and in the far SOL where, in certain cases, the formation of a shonlder in the density profile is 
observed. The experiments that we performed showed that the far SOL flattening can occur 
even in the absence of divertor detachment or wall recycling, although both mechanisms are 
likely to play a role in particular cases. It is interesting that even 2D turbulence models 
without divertor or neutral particle physics capture the prohle broadening 
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30l |. although 


they also suggest that a proper 3D treatment in realistic geometry is required to unravel the 


details behind it. W 
only recently 
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rile the 3D dynamics of the hlaments were addressed by the community 


39| . the hrst results already anticipate a rich phenomenology in which 


both local and global features of the plasma hlaments can affect their motion in the SOL. 
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In particular, crucial effects could include the presence and structure of the X-point; the 
parallel prohle of the temperature (or collisonality) in the hlament; the electrical connection 
of the hlament to the Debye sheath at the target through the background plasma or its 
absence; the parallel and perpendicular shape and dimension of the hlament; the ionization 
in the hlament caused by its higher temperature with respect to the background. While 
the divertor regime (its collisonality, the occurrence of detachment) certainly ahects some 
of the items above, it does not provide a complete picture. Similarly, wall recycling changes 
the particle source, which might be important for machines with a narrow gap between 
separatrix and wall, but does not have the generality to capture the phenomena observed in 
MAST. 

It is interesting to observe that the shoulder formation is not observed in our high current 
and low density low magnetic held discharges, which suggests that the safety factor could 
play a role in this phenomenon. On the other hand, we do not rule out the possibility that 
the broadening can occur at larger distances from the separatrix in these latter cases (i.e. 
the far SOL is farther away from the separatrix). In our experiments we did not monitor 
the SOL at distances beyond ~ 5cm as this would imply analysing reciprocating probe data 
corresponding to a diherent plasma regime (since the probe plunges at a certain velocity in 
the plasma, the far SOL corresponds to earlier times in the discharge). From a practical 
point of view, however, this makes little difference, since the density prohle would anyway 
reach low levels before battening, thus limiting the plasma wall interactions. 

From an experimental point of view, pinpointing what exactly causes the battening on 
the near and far SOL proved extremely difbcult. This was due to the limited database 
and the uncertainty in the measurements, notably the signibcant errorbars on the position 
of the separatrix (~ 1cm) and on the edge temperature (relative error < 65%). Both 
these quantities abect the determination of the collisonality ~ which is a scaling 

parameter often used to interpret and extrapolate the SOL particle exhaust. While more 
precise measurements of the upstream temperature are desirable and to a certain extent 
already possible using other diagnostics, the quadratic dependence of the collisonality on 
this quantity makes it problematic anyway. In addition, another degree of arbitrariness is 
introduced by the singular behaviour of the safety factor on the last closed bux surface, 
which allows a continuum of connection lengths (and therefore collisionalities) depending 
on the radial distance from the separatrix where this quantity is estimated. In order to 
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resolve this conundrum, only the combination of better measurements with hrst principle 
understanding of the SOL mechanisms is a viable solution. 

Differently from the particle exhaust, the energy exhaust displays a remarkable insensi¬ 
tivity to the thermodynamic state of the main plasma, at least at the target. This is in 
agreement with experimental multi-machine scaling laws [l^, which relate the parallel heat 
flux decay length mainly to the magnetic structure of the plasma. As a matter of fact, our 
target data are well represented by the L-mode Scarabosio scaling. The weak dependence 
of the target decay length on the upstream density and temperature is further conhrmed 
by its response to sawtooth oscillations. We consistently observed a self-similar prohle in 
the divertor IR data during the sawtooth cycle, apart from a short period around the crash 
when the target prohle has to dynamically respond to the higher heat hux received. Our 
analysis cannot clarify if the absence of thermodynamic effects on the target Xq implies 
that magnetic geometry dominates over plasma physics or that there is an almost complete 
cancellation of the latter. The upstream heat hux decay length seems to be more sensitive, 
albeit only moderately, to the local density and temperature and are reproduced within error 


bars by theoretical scaling laws based on upstream turbulence 


30| . The inner target heat 


hux deposition features a Gaussian prohle rather than the convolution between a Gaussian 
and an exponential used in the outer divertor. This suggests a quite diherent in/out exhaust 
mechanism, at least in the double null conhgurations examined. 

The results above suggest that the particle and energy exhaust occur through diherent, 
albeit probably related, mechanisms. This is further conhrmed by the fact that the midplane 
projection of the divertor data, based on perpendicular dihusion models, gives a good match 
with the measured prohles of the heat hux but not of the saturation current (i.e. the 
density). A possible explanation for this discrepancy relies on the diherent parallel transport 
paradigms followed by the particle and energy channels j^|. The heat is efficiently exhausted 
along held lines by the electrons, which are highly conductive in the parallel direction. This 
means that, whatever the mechanism that moves the energy across held lines (turbulent 
eddies, hlaments, collisions), it has little time to act before most of the heat is drained 
at the divertor. Gonversely, the particle exhaust in the parallel direction is advective and 
much slower, so that the characteristic length scale for the SOL density, in the range of 
collisionality considered, is larger than that of the heat hux. In conclusion, as the steep 
heat hux prohle samples mainly the dihusive region close to the separatrix, the projection 
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technique works well, while it does not for the particle exhaust, characterised by non-local 
non-diffusive transport related to the filaments. 

Finally, the limited number of discharges performed did not allow us to extract reliable 
scaling laws from the data collected. This explains the mostly qualitative nature of our 
discussion and advocates for multi-machine analysis in order to provide quantitative relations 
between the relevant parameters of the problem. 


VIII. ACKNOWLEDGEMENTS 

F.M. acknowledges useful discussions with Prof. B. Lipschultz, Dr. S. Elmore and Mr. 
N. Walkden. We aslo thank Dr. K. Lawson for carefully reading the manuscript. This work 
has been carried out within the framework of the EUROfusion Consortium and has received 
funding from the Euratom research and training programme 2014-2018 under grant agree¬ 
ment No 633053 and from the RCUK Energy Programme [grant number EP/I501045]. To 
obtain further information on the data and models underlying this paper please contact Pub- 
hcationsManager@ccfe.ac.uk. The views and opinions expressed herein do not necessarily 
reflect those of the European Commission. 


[1] A.Loarte, B. Lipschultz, A.S. Kukushkin et ai, Nucl. Fusion 47, S203 (2007). 

[2] P. C. Stangeby, The Plasma Boundary of Magnetic Fusion Devices, (CRC Press, 2000). 

[3] G.Y. Antar, S. 1. Krasheninnikov, P. Devynck, R. P. Doerner, E. M. Hollmann, J. A. Boedo, 
S. C. Luckhardt, and R.W. Conn, Phys. Rev. Lett. 87, 065001-1 (2001). 

[4] G. Y. Antar, G. Counsell, Y. Yu, B. Labombard, and P. Devynck, Phys. Plasmas 10, 419 
(2003). 

[5] B. LaBombard, M.V. Umansky, R.L. Boivin, J.A. Goetz, J. Hughes, B. Lipschultz, D. Mosses- 
sian, C.S. Pitcher, J.L. Terry, Alcator Group, Nucl. Fusion 40, 2041 (2000). 

[6] B. LaBombard, R. L. Boivin, M. Greenwald, J. Hughes, B. Lipschultz, D. Mossessian, C. S. 
Pitcher, J. L. Terry, S. J. Zweben, Phys. Plasmas 8, 2107 (2001). 

[7] B. Lipschultz, B. LaBombard, C.S. Pitcher and R. Boivin, Plasma Phys. Control. Fusion 44, 
733 (2002). 

28 



[8] D.G. Whyte, B.L. Lipschultz, P.C. Stangeby, J. Boedo, D.L. Rudakov, J.G. Watkins and W.P. 
West, Plasma Phys. Gontrol. Fusion 47, 1579 (2005). 

[9] O.E. Garcia, J. Horacek, R.A. Pitts, A.H. Nielsen, W. Fnndamenski, V. Nanlin and J. Junl 
Rasmnssen, Nucl. Fnsion 47, 667 (2007). 

[10] D. Garralero, G. Birkenmeier, H.W. Muller, P. Manz, P. deMarne, S.H. Muller, F. Reimold, 
U. Stroth, M. Wischmeier, E. Wolfrum, Nucl. Eusion 54, 123005 (2014). 

[11] T. Eich, B. Sieglin, A. Scarabosio et ah, Phys. Rev. Lett. 107, 215001 (2011). 

[12] M. A. Makowski, D. Elder, T.K. Gray, B. LaBombard, G.J. Lasnier, A. W. Leonard, R. 
Maingi, T.H. Osborne, P.C. Stangeby, J.L. Terry, and J. Watkins, Phys. Plasmas 19, 056122 
( 2012 ). 

[13] T. Eich, A.W. Leonard, R.A. Pitts, W. Eundamenski, R.J. Goldston, T.K. Gray, A. Herrmann, 
A. Kirk, A. Kallenbach, O. Kardaun et al., Nucl. Eusion 53, 093031 (2013). 

[14] A. Scarabosio, T. Eich, A. Herrmann and B. Sieglin, J. Nucl. Mat. 438, S426 (2013). 

[15] E. Militello and W. Eundamenski, Plasma Phys. Control. Eusion 53, 095002 (2011). 

[16] A. Sykes, R.J. Akers, L.C. Appel, E.R. Arends, P.G. Carolan, N.J. Conway, G.E. Counsell, 
G. Cunningham, A. Dnestrovskij, Yu.N. Dnestrovskij et al, Nucl. Eusion 41, 1423 (2001). 

[17] J.R. Harrison, S.W. Lisgo, K.G. Gibson, P. Tamain, J. Dowling, J. Nucl. Mat. 415, S379 

( 2011 ). 

[18] G. De Temmerman, E. Delchambre, J. Dowling, A. Kirk, S. Lisgo and P. Tamain, Plasma 
Phys. Control. Eusion 52, 095005 (2010). 

[19] R. Scannell, M. J. Walsh, M. R. Dunstan, J. Eigueiredo, G. Naylor, T. OGorman, S. Shibaev, 
K. J. Gibson, and H. Wilson, Review of Scientific Instruments 81, 10D520 (2010). 

[20] Y. Yang and G. Counsell, Journal of Nuclear Materials 313316, 734 (2003). 

[21] L.L. Lao, H. St. John, R.D. Stambaugh, A.G. Kellman and W. Pfeiffer, Nuclear Eusion 25, 
1421 (1985). 

[22] S. Elmore, S.Y. Allan, A. Kirk, G. Eishpool, J. Harrison, P. Tamain, M. Kocan, R. Gaffka, R. 
Stephen, J.W. Bradley, Plasma Phys. Control. Eusion 54, 065001 (2012). 

[23] G.M. McCracken, M.E. Stamp, R.D. Monk, A.G. Meigs, J. Lingertat, R. Prentice, A. Starling, 
R.J. Smith, A. Tabasso, Nucl. Eusion, Vol. 38, 619 (1998). 

[24] J. Huang, S. Lisgo, G. Maddison, Plasma Phys. Control. Eusion 52, 075012 (2010). 


29 


[25] K. McCormick, G. Kyriakakis, J. Neuhauser, E. Kakoulidis, J. Schweinzer and N. Tsois, J. 
Nucl. Mat. 196, 264 (1992) 

[26] N. Asakura, Y. Koide, K. Itami, N. Hosogane, K. Shimizu, S. Tsuji-Iio, S. Sakurai, A. Sakasai, 
J. Nucl. Mat. 241, 559 (1997). 

[27] B. LaBombard, J.W. Hughes, D. Mossessian, M. Greenwald, B. Lipschultz, J.L. Terry, Nucl. 
Fusion 45, 1658 (2005). 

[28] B. Lipschultz, D. Whyte and B. LaBombard, Plasma Phys. Control. Fusion 47, 1559 (2005). 

[29] F. Militello, W. Fundamenski, V. Naulin and A.H. Nielsen, Plasma Phys. Control. Fusion 54, 
095011 (2012). 

[30] F. Militello, V. Naulin and A.H. Nielsen, Plasma Phys. Control. Fusion 55, 074010 (2013). 

[31] J.T. Omotani and B.D. Dudson, Plasma Phys. Control. Fusion 55, 055009 (2013). 

[32] R.J. Goldston, Nucl. Fusion 52, 013009 (2012). 

[33] F. Wagner, Nucl. Fusion 25, 525 (1985). 

[34] F. Militello et al, in preparation (2015). 

[35] J.R. Angus, M.V. Umansky, and S.I. Krasheninnikov, Phys. Rev. Lett. 108, 215002 (2012). 

[36] J.R. Angus, M.V. Umansky, and S.I. Krasheninnikov, Phys. Plasmas 19, 082312 (2012). 

[37] F.D. Halpern, A. Cardellini, P. Ricci, S. Jolliet, J. Loizu and A. Mosetto, Physi. Plasmas 21, 
022305 (2014). 

[38] L. Easy, F. Militello, J. Omotani, B. Dudson, E. Havhckova, P. Tamain, V. Naulin, and A. 
H. Nielsen, Phys. Plasmas 21, 122515 (2014). 

[39] J. Omotani, F. Militello and L. Easy, in preparation (2015). 

[40] W. Fundamenski, O.E. Garcia, V. Naulin, R.A. Pitts, A.H. Nielsen, J.J. Rasmussen, J. Ho- 
racek, J.P. Graves and JET EFDA contributors, Nucl. Fusion 47, 417 (2007). 


30 


10'° D2/m3 


Shot: - 28996 - 28998 - 29003 - 29007 - 29018 - 29023 - 29027 

SMTH(XIM_DA/HM10/R, 0.01 SEC) 



4.0 



0.20 0.22 0.24 0.26 0.28 0.30 0.32 


FMILITEL Wed Dec 3 16:07:30 2014 


Time (Sec) 


XPAD6 (Version 3E - Idem) 



















